Accelerating equilibrium isotope effect calculations: I. Stochastic thermodynamic integration with respect to mass
نویسندگان
چکیده
Accurate path integral Monte Carlo or molecular dynamics calculations of isotope effects have until recently been expensive because of the necessity to reduce three types of errors present in such calculations: statistical errors due to sampling, path integral discretization errors, and thermodynamic integration errors. While the statistical errors can be reduced with virial estimators and path integral discretization errors with high-order factorization of the Boltzmann operator, here we propose a method for accelerating isotope effect calculations by eliminating the integration error. We show that the integration error can be removed entirely by changing particle masses stochastically during the calculation and by using a piecewise linear umbrella biasing potential. Moreover, we demonstrate numerically that this approach does not increase the statistical error. The resulting acceleration of isotope effect calculations is demonstrated on a model harmonic system and on deuterated species of methane.
منابع مشابه
Direct path integral estimators for isotope fractionation ratios.
Fractionation of isotopes among distinct molecules or phases is a quantum effect which is often exploited to obtain insights on reaction mechanisms, biochemical, geochemical, and atmospheric phenomena. Accurate evaluation of isotope ratios in atomistic simulations is challenging, because one needs to perform a thermodynamic integration with respect to the isotope mass, along with time-consuming...
متن کاملPath integral evaluation of equilibrium isotope effects.
A general and rigorous methodology to compute the quantum equilibrium isotope effect is described. Unlike standard approaches, ours does not assume separability of rotational and vibrational motions and does not make the harmonic approximation for vibrations or rigid rotor approximation for the rotations. In particular, zero point energy and anharmonicity effects are described correctly quantum...
متن کاملHydration in solution is critical for stable oxygen isotope fractionation between carbonate ion and water
Various isotope studies require accurate fractionation factors (a’s) between different chemical compounds in thermodynamic equilibrium. Although numerous isotope systems involve aqueous solutions, the conventional theory is formulated for the gas-phase and predicts incorrect a’s for many compounds dissolved in water. Here I show that quantum-chemistry calculations, which take into account solut...
متن کاملModeling Basel Regulatory in DSGE with Emphasis on Adequacy Regulatory
In this paper Basel regulation is modeled in Dynamic Stochastic General Equilibrium (DSGE) framework. For this purpose, using data from 1981-2017 for Iran, capital adequacy as an importance regulation is modeled. Results show Basel regulation has procyclical effect. According to the results of the model and according to the realities of economy and banking system of Iran, in recession, lending ...
متن کاملInvestigating the Effects of Mass Transfer and Mixture Non-Ideality on Multiphase Flow Hydrodynamics Using CFD Methods
A numerical framework has been proposed to model the interacting effects of mixture non-ideality and mass transfer on hydrodynamics of a multiphase system using CFD methods.Mass transfer during condensation and vaporization is modeled by chemical potential at the liquid-vapor interface. Species mass transfers are related to the diffusion at the interface which in turn is related to the conc...
متن کامل